Granulated superconductors: from the nonlinear a model to the Bose-Hubbard 

description 
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We modify a nonlinear a model (NLaM) for the description of a granulated disordered system 
in the presence of both the Coulomb repulsion and the Cooper pairing. We show that under certain 
controlled approximations this model is reduced to the Bose-Hubbard (or "dirty-boson" ) model with 
renormalized coupling constants. We obtain a more general effective action (which is still simpler 
than the full NLctM action) which can be applied in the region of parameters where the reduction 
to the Bose-Hubbard model is not justified. This action may lead to a different picture of the 
superconductor-insulator transition in 2D systems. 
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A wide variety of experimental data on the super- 
conductor-insulator (SI) transition in two-dimensional 
structures continues to attract acute theoretical in- 
terest. The transition can be tuned by either disorder 
(changing with the thickness of a superconducting film) 
or magnetic field, thus being one of the most intensely 
studied examples of quantum phase transitions [|| . How- 
ever, recent experiments have challenged the very 
existence of the SI transition, leaving open a possibility 
that a dramatic drop in resistance is due to the existence 
of a crossover to a new metallic phase with resistance 
much lower than that in the normal state and to a subse- 
quent metal-superconducting transition This situa- 
tion requires a reassessment of theoretical approaches to 
the problem of dirty superconductors. 

One of the ways to understand the problem of the 
SI transition is based on the so-called Bose-Hubbard 
(or "dirty-boson" ) models Q where the superconducting 
phase is due to the bose-condensation of charge-2e bosons 
(preformed Cooper pairs) with localized vortices while 
the insulating phase is due to the bose-condensation of 
vortices with localized Cooper pairs. Another approach 
which captures the basic physics of granular supercon- 
ductors is based on dissipative modelsj^ of resistively 
shunted charged Josephson arrays [p[-|ll [, with the em- 
phasis on the role of dissipation and Coulomb interac- 
tion. In both group of models [^-|Tl|], the transition 
is driven by fluctuations of the phase of the order pa- 
rameter. An alternative approach is based on a micro- 
scopic description of homogeneous systems that incorpo- 
rates both the attractive (in the Cooper channel) and 
repulsive electron-electron interaction in the presence of 
disorder into an effective field theory, the nonlinear a 
model (NLaM) . The SI transition in these mod- 

els is driven by fluctuations of the amplitude rather than 
the phase of the order parameter, and the Cooper pairing 
is suppressed by the repulsive interaction on the insulat- 
ing side of the transition. 

The purpose of this paper is to show that both 
the Bose-Hubbard model and the dissipative models 



" |11[ can be in fact derived from the same microscopic 
action of a NLctM modified for the description of granu- 
lar superconductors. The two models correspond to cer- 
tain simplifications made within the NLctM. The latter is 
more general and allows one to go beyond different limi- 
tations inevitable in the derivation of the Bose-Hubbard 
and dissipative models. Note that the dissipative action 
1^ of Ambegaokar, Eckern, and Schon (AES) has been 
widely used p^ ] in a simplified form in the context of a 
normal tunnel junction. This variant of the AES action 
has been very recently derived ijl^ from the NLaM de- 
scribing electrons with the repulsive interaction moving 
in the presence of disorder. Here we will derive both the 
full AES action for Josephson junctions and the Bose- 
Hubbard model from the NLcrM that includes both the 
attraction in the Cooper channel and the Coulomb repul- 
sion. We shall use a new variant ||l^ of the NLcrM which, 
in our opinion, considerably simplifies the calculations. 
Naturally, one can use (after a straightforward modifi- 
cation for a granular system) any version of the NLcrM 
that includes the Cooper pairing and the Coulomb in- 
teraction, either the original Finkelstein model ||T^ or a 
more recent model in Keldysh technique. 

Our starting point is the standard microscopic Hamil- 
tonian which includes a (5-correlated Gaussian random 
potential, the Coulomb interaction and the BCS attrac- 
tion. We consider a coarse-grained version of this Hamil- 
tonian which corresponds to a granular superconductor. 
This will allow us to separate scales of fluctuations of 
the amplitude and the phase of the superconducting or- 
der parameter A. Neglecting (at some later stage) the 
amplitude fluctuations and making some further simpli- 
flcations will lead eventually to the models |^,^ governed 
only by the phase fluctuations. 

Let us stress again that the NLcrM used here can be 
applied to the simultaneous description of both the am- 
plitude and phase fluctuations of A that can be quite im- 
portant in the relation to recent experiments [0. More- 
over, even when focusing on the phase fluctuations only, 
the effective functional is essentially generalized by dis- 
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order (affecting intra-granular electron motion and thus 
leading to a different model of the phase fluctuations) 
and can in principle lead to a different picture of the 
transition. 

The derivation of the NLcM from the Hamiltonian 
described above follows the standard steps First 
one averages the replicated imaginary-time fermionic ac- 
tion over the random potential. Then one employs the 
Hubbard- Stratonovich transformation to decouple three 
quartic (in the electron field) terms, corresponding to the 
disorder, the Coulomb repulsion and the BCS attraction. 
Finally, integrating out the fermionic fields, one arrives 
at the effective action in terms of three bosonic fields: a 
matrix field a which decouples the disorder-induced "in- 
teraction", $ which decouples the Coulomb repulsion, 
and A which decouples the BCS attraction: 
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Here the operator e equals if^d-r in imaginary time repre- 
sentation and becomes the diagonal matrix of fermionic 
Matsubara frequencies in frequency representation, ^ is 
the operator of the intra-grain kinetic energy (counted 
from the chemical potential), and t is the tunneling am- 
plitude matrix (i.e. the inter-grain kinetic energy). All 
the bosonic fields are defined in the space which is con- 
venient to think of as a direct product of the N x N 
replica sector, 2x2 spin sector, 2x2 'time-reversal' sec- 
tor (introduced for a correct decoupling in the Cooper 
channel for both disorder- induced and BCS interaction), 
and of the m x m grain sector. The symbol Tr refers 
both to a summation over all these matrix indices and 
to an integration over position r and the imaginary time 
T. The matrix a is diagonal in grain indices (called later 
and possesses standard symmetries in all the other 
sectors ||T^ . The pairing field A is diagonal in the replica 
and grain indices and in x = (r, t), and has the following 
structure |l8| in time-reversal and spin sectors: 



A{x) = |A(x)|e^x(:'=)?3^«P 
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where Tq and are Pauli matrices in the time-reversal 
and spin sectors, respectively. The Coulomb field ^{x) 
is proportional to the unit matrix in all the matrix sec- 
tors. Finally, in Eq. (|l|) Ao and t^i are the BCS coupling 
constant and the elastic mean free time, respectively, and 
U = U{r — r') is the Coulomb interaction. 

The principal simplification for granular systems is 
that all the fields are spatially homogeneous inside each 
grain when the grains are zero-dimensional, i.e. their 
sizes L ^ ^, Lt {£, and Lt are the superconducting 
and thermal coherence lengths) which is equivalent to 
|A|,T ^ 1/Tcrg- Then the Coulomb interaction reduces 
to the capacitance matrix, — > Cij/e"^, and the tun- 
neling matrix t — {tij} depends only on grain indices. 



Now we follow the procedure of derivation of the 
NLctM for dirty superconductors First, we look for 
a saddle point of (||) with respect to a separately for each 
grain. It can be parameterized as as. p. = S^AS, where S 
is a certain matrix which also diagonalizes e-|-<I>-|-A, 



Ai — SjXiSi , 



(3) 



where all the matrices are diagonal in grain indices. Then 
the entire saddle-point manifold is parameterized as 
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where in the Matsubara representation A = diagjsgne}, 
and matrix U defines the standard coset space po|] . Here 
Q is a degenerate solution to the saddle-point equation 
for the action (|l|) when e = 0, A and <& all vanish. The 
parameterization (^ 'aligns' the field a so that A and <i> 
are taken into account in the zeroth approximation. 

Now we perform a similarity transformation with ma- 
trices S and under Trln in Eq. (^. As all the fields 
are spatially homogeneous inside each grain, S commutes 
with the operator ^. Then one only needs to expand the 
Tr In to the first nonvanishing orders in tij and , this 
expansion being justified when |A|,r ^ 1/tci ^ Ef- 
Thus one arrives at the following effective action: 
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^Tr XiQi ^ ^ Tr Q^SijQjSji , 



(5) 



where Sij = SiSj, all the fields depend on r, Tr refers 
to all indices except those numerating grains, 6i is mean 
level spacing in the i-th grain, and the tunneling conduc- 
tance is defined by gfj = 2TT'^\tij\'^ /6i6j (which is nonzero 
only for neighboring grains). Both S and A should be 
found from the diagonalization procedure in Eq. (H) . 
The next step is to represent Si as 



S, = K-e"*^'(^)^=' 



(6) 



This is similar to the gauge transformation suggested in 
Refs. and used in Ref. to gauge out the Coulomb 
field. However, one cannot gauge out two independent 
fields, A and Substituting the transformation (||) into 
the diagonalization condition (^, we reduce it to 



A° = V+X,V, , 



(7) 



where A° is the field (g) taken at x = and the field $ 
is given in the t representation by $j = <i>i ~ ^drXi ■ 

Both <i> and |A| are massive fields whose fluctuations 
are strongly suppressed. It is straightforward to show 
that the fluctuations of $ are of order S which is much 
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smaller than both T and | A|. Therefore, in the mean field 
approximation in $ this field can be neglected, $ = 0. 
This condition is nothing' more than the Josephson rela- 
tion in imaginary time Q . This locks the fluctuations of 
the Coulomb field <f> with the phase-fluctuations of the 
pairing field A, 



(8) 



thus reducing the action (|^) to one depending only on 
the fields Q and A. 

The mean field approximation in |Ai| is valid for 
I Ai I » (5 and and reduces to the standard self-consistency 
equation which formally follows from the variation of the 
action (^) with respect to |A|. In this approximation one 
finds I Ail to be independent of i and r. Thus the first 
term in Eq. becomes a trivial constant, so that the 
action depends on Q and % only: 
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Note that the field x iii fl^is action obeys the standard 
boundary condition x(t + /?) = xif) mod 2-k. Thus when 
calculating the partition function with this action, one 
should take into account different topological sectors cor- 
responding to different winding numbers in x- 

Now the diagonalization conditions ^ are the same 
for each grain, and reduce to those solved in Ref. pi]: 



Vee' = cos Y -|- ®T2 siu y SgUe 5^^^^, 
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A = diag-v/ £2 + |Ap sgn e , cos 9^ = — , . 
Then Sij in Eq. (^) can be expressed in terms of V as 



XiJ =Xr- Xj ■ 



(11) 



Finally note that large- |e| contributions to the action (Q) 
are strongly suppressed, while for |e| <C |A| one has 
A = |A|A which suppresses fluctuations of Q in each 

fdn imposing Q = A. Then, all matrices in the action 
are diagonal in the replica indices so that these in- 
dices become redundant. The diagonalization procedure 
(|l^)-(^l]) has resolved explicitly the matrix dependence 
on the time-reversal and spin indices. This reduces the 
action to that depending only on one scalar bosonic field, 
the phase x of the order parameter, whose arguments are 
the imaginary time r, and the position index i, i.e. the 
grain number. Indeed, the second term in Eq. ^ reduces 
to a trivial constant; evaluating the tunneling term with 
the help of Eqs. (|lO|) and ([l|), we obtain 
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where Xjj(t) = Xiir) - Xi(-r), 



Xij (t) ± Xij (t') 



and the normal and anomalous Green's functions gn^a 
(integrated over all momenta) are given by 



5«(t) = T 



""^^^ :,5.(r) = TV^I^. (13) 



Ve2 + |A|2 



This action coincides with that derived in Refs. 
Further simplifications are possible in two limiting cases 
First, in the normal case (A = 0), one has in Eq. (13) 



5a = 0, 



glir) 



rp2 



sin ttTt 
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Then the field x should be substituted, according to 
Eq. (^), by 2 /^dr'^jr'). This hmiting case corresponds 
to using the action (|lj) , (|lj) in the context of a normal 
tunnel junction p6[ |. This is precisely the action which 
has been recently derived from the NLctM in Ref. [ pT| ; 
the functional ^ in the limit A = is equivalent to 
the cr model of Ref. . Including the disorder- induced 
fluctuations (i.e. going beyond the Q = A approxima- 
tion) allows one to obtain [jl^ a correct low-T limit for 
the phase correlation function missing in the action ([l^). 

The action (|^) is more general than that considered 
in Ref. [0: although under the mode locking condition 
(^) it depends only on the fields x a-nd Q, the matrix 
Sij, Eqs. ( |lO| ) and ([ll|), reduces to a simple U{1) gauge 
transformation as in in Ref. only in the hmit A = 0. 

The second limiting case, T <^ |A|, is just the limit 
relevant in the context of the SI transition in granular 
superconductors. For T = 0, the summation in Eq. ( p^ ) 
can be substituted by integration which yields 

g„(r) = ^ifi(|A|r), g^{r) ^ Ko{\A\t) 

This is also a good approximation for a low-temperature 
case; substituting this into Eq. ( p^ ) gives the action for 
the dissipative model [^,0. Note that for |t — t'| ^ 



A|^^, the main contribution in the tunneling action (|1^ 
is given by the normal term with the corresponding ker- 
nel proportional to \t ~ t'\^'^ . The Fourier transform of 
this would give a term of the Caldeira-Leggett type 0| 
proportional to 

The tunneling action ( |l^ ) is non-local in r. As has 
been noted in Ref. for the case of one tunnel junction, 
for sufficiently large capacitance the phase Xij changes 
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slowly in comparison with |A|~"'^, and in the adiabatic 
approximation x(''"') is changed by x(t) + (r'— t)9t-x(t). 
Making such an expansion, one obtains from Eq. (^2|) the 
following local action: 



" 

where Xi = drXi ^-nd 
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gjj 3 + cos Xij 
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_i ^ _ fft-,- 3 + cos Xij 
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considerably different from the NLcrM action for homo- 
geneous systems. We hope that using this action may 
eventually lead to a different phase diagram for granu- 
lated superconductors. 
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If all the self-capacitances are equal to C with Ec oc e^/C 
being the charging energy, and all gjj — , then uu = U 
has the meaning of the renormalized charging energy. Ig- 
noring a weak dependence of u on cos Xij in the above 
relations, one obtains the renormalized charging energy: 



U ■ 



Er 
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(16) 



Here the coefficient jj^ depends on the number of next 
neighbors for each grain, etc. A similar renormaliza- 
tion takes place for the next-neighbor off-diagonal energy 
Uij . Now one can see that on the face of it the adiabatic 
approximation employed to obtain Eq. (Rq) is valid for 
U <^ |A|. However, in the region ^\A\/Ec, where 
the charging energy (|lq) is strongly renormalized, the 
instanton-like solutions [22[| may be important. This may 
further reduce the region of applicability for the local in 
T action (15). 

Finally, by introducing the operator n canonically con- 
jugate to the phase x, one finds the Hamiltonian that 
corresponds to the action ([15|): 



H = 



HI 



2 ^ij^i^j 



^\glj cos (xi - X]] 



(17) 



This is just the Hamiltonian of the Bose-Hubbard model 
which was first microscopically derived by Efetov | |2^ ] 
in the context of granulated superconductors. 

To conclude, we have shown that for a granular sys- 
tem with zero-dimensional grains the NLcrM reduces to 
the AES model (|l^). When the charging energy is 
much smaller than |A|, the action (n2|) further reduces 
to the Bose-Hubbard model, Eq. (|17|), which is widely 
used for the description of the superconductor-insulator 
transition However, the above estimations show 

that this reduction is parametrically justified only for 
the region Ec <C |A| where the transition happens at 
g"^ -C Ec/\A\ <C 1 which corresponds to a strongly gran- 
ulated system. Note finally that the most general (in 
the present context) action (^ describes both amplitude 
and phase fluctuations of the order parameter, being still 
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